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ABSTRACT 



With the 100-m telescope at Effelsberg, 19 ammonia (NH3) maser lines have been detected toward the prominent 
massive star forming region W51-IRS2. Eleven of these inversion lines, the (J \K) = (6,2), (5,3), (7,4), (8,5), (7,6), 
(7,7), (9,7), (10,7), (9,9), (10,9), and (12,12) transitions, are classified as masers for the first time in outer space. All 
detected masers are related to highly excited inversion doublets. The (5,4) maser originates from an inversion doublet 
~340K above the ground state, while the (12,12) transition, at -1450K, is the most highly excited NH3 maser line 
so far known. Strong variability is seen not only in ortho- but also in para-NH3 transitions. Bright narrow emission 
features are observed, for the first time, in (mostly) ortho-ammonia transitions, at VLsr ~ 45kms , well separated 
from the quasi-thermal emission near 60kms _1 . These features were absent ~25 years ago and show a velocity drift of 
about +0.2kms _1 yr _1 . The component is likely related to the SiO maser source in W51-IRS2 and a possible scenario 
explaining the velocity drift is outlined. The 57kms _1 component of the (9,6) maser line is found to be strongly linearly 
polarised. Maser emission in the (J,K) to (J+1,K) inversion doublets is strictly forbidden by selection rules for electric 
dipole transitions in the ground vibrational state. However, such pairs (and even triplets with (J+2,K)) are common 
toward W51-IRS2. Similarities in line widths and velocities indicate that such groups of maser lines arise from the same 
regions, which can be explained by pumping through vibrational excitation. The large number of NH3 maser lines in 
W51-IRS2 is most likely related to the exceptionally high kinetic temperature and NH3 column density of this young 
massive star forming region. 
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1. Introduction 



ISM: clouds - ISM: individual objects: W51 - ISM: HII regions - ISM: molecules - Radio lines: 



Ammonia (NH3) provides the unique opportunity to trace 
molecular cloud excitation up to temperatures of — 2000K 
by observing its characteristic inversion transitions within 
a very limited frequency interval (20 - 35 GHz; e.g. Ho & 
Townes 1983; Wilson et al. 2006, 2008). The inversion dou- 
blets arise from oscillations of the nitrogen nucleus through 
the plane of the three hydrogen nuclei. The frequencies of 
the lines connecting the two states of an inversion-doublet 
depend on the total angular momentum J and its projec- 
tion on the molecular axis, K, with K = 0,3,6,9... belonging 
to ortho-NH3 and K = 1,2,4,5,7... representing para-NH3. 
Thus in the vibrational ground state, inversion doublets are 
characterized by the values J and K. The parities of the 
split levels differ to allow for the dipole transitions (J,K) — > 
(J,K) (hereafter (J,K)), which are subject of this paper. 

Dozens of inversion lines can be detected, provided ki- 
netic temperatures and ammonia column densities are high 
enough. These conditions prevail in "hot cores", dense 
molecular clumps near sites of very recent massive star 
formation (e.g., Mauersberger et al. 1986a, 1988b; Henkel 
et al 1987b; Hermsen et al. 1988; Cesaroni et al. 1992; 
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Hiittemeister et al. 1993, 1995; Wilson et al. 1993, 2000; 
Zhang & Ho 1997; Goddi et al. 2011). The NH 3 abundances 
are believed to be caused by dust grain mantle evapora- 
tion (e.g., Henkel et al. 1987a; Walmsley et al. 1987; Brown 
et al. 1988). The warm dense clumps are characterized by 
temperatures T kin > 100 K, AT(NH 3 ) = AT(NH 3 )/./V(H 2 ) ~ 
2Q-5...-6 an( j source averaged ammonia column densities in 
excess of 10 18 cm~ 2 . 

Most but not all the ammonia inversion lines are ther- 
mally excited. (J,K) = (3,3) maser emission was first de- 
tected by Wilson et al. (1982) toward the star forming re- 
gion W33. To date, several NH3 maser lines have been ob- 
served (e.g., Guilloteau et al. 1983; Madden et al. 1986; 
Mauersberger et al. 1988b; Hofner et al. 1994; Kraemer & 
Jackson 1995; Beuther et al. 2007), even including the rare 
isotopologue 15 NH 3 (Mauersberger et al. 1986b; Schilke et 
al. 1991). 

A particularly outstanding hot core is associated with 
the star forming region W51-IRS2, which includes W51- 
North, W51d, W51dl, and W51d2 (Zapata et al. 2009; their 
Fig. 1). At a distance of 6.1±1.3kpc (Imai et al. 2002) or 
5.lt^kpc (Xu et al. 2009), lpc subtends -40". With a 
rotational temperature T rot — 300 K, a source averaged col- 
umn density of A^(NH 3 ) — 10 19 cm -2 , and a virial density 
of n(H 2 ) ~ 10 7 7 cm- 3 , W51-IRS2 is the source with the 
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largest number of detected masing ammonia lines (Madden 
et al. 1986; Mauersberger et al. 1987; Wilson & Henkel 
1988). Here we focus on such non-thermally excited lines 
and report results from monitoring observations spanning 
a time interval of almost two decades and leading to the 
detection of a dozen new molecular maser lines. 



2. Observations 

The data were taken in November and December 1995, 
October 2008, August 2011, and in April 2012 with the 
100-m Effelsberg telescop43 near Bonn/Germany. At the 
line frequencies observed (18.5 - 31.4 GHz), the beam size 
is 49" - 29". The main beam brightness temperature scale 
was established by continuum cross scans toward 3C 286 
and NGC 7027 (flux densities were adopted from Ott et al. 
(1994), also accounting for a 0.5% yr _1 secular decrease in 
the case of NGC 7027), while 1923+210 and 2145+06 were 
used as pointing sources. The pointing accuracy was bet- 
ter than 10". For maser sources much more compact than 
the beam size, line shapes are not affected by pointing er- 
rors. In 1995 and 1996, at frequencies above 26 GHz, a sin- 
gle channel primary focus HEMT (Hot Electron Mobility 
Transistor) receiver was employed with a system equiva- 
lent temperature on a Jansky scale of ~ 85 Jy, including 
sky and ground radiation. Later, measurements at K-band 
(18-26 GHz) were carried out with a dual channel (includ- 
ing both orthogonal linear polarizations) cooled primary 
focus HEMT receiver with a T sys equivalent of ~65 Jy per 
channel. 

For the early observations, the backend was a three level 
autocorrelator with 1024 spectral channels. Applied band- 
widths ranged from 3.125 to 25 MHz. In 2008, and 2009, 
a Fast Fourier Transform Spectrometer (FFTS) was em- 
ployed with 16384 channels for each polarization. Applied 
bandwidths were 50 and 500 MHz. For the most recent data 
another FFTS was used, this time with a bandwidth of 
2 GHz and 32768 channels per polarization. This allowed 
us to cover the entire 18-26 GHz range with a few fre- 
quency setups. Following weather conditions in a flexible 
way (measurements are much less weather dependent at 
18-20 GHz than at higher frequencies), the frequency setup 
was changed quite often during the most recent observa- 
tions with the 2 GHz backend. These covered parts of five 
consecutive days. Therefore, averaging data from several of 
these days, April 4-8, 2012, is mostly given as the epoch 
for these spectra. 

All data were taken in the position switching mode, em- 
ploying scans of 6 minutes total duration and offset po- 
sitions 10' east or west. The calibration uncertainties are 
estimated to be ±15% (lc). 



3. Results 

3.1. Are there lines caused by maser action? 

Figures HH7] present spectra taken during a time interval 
of almost two decades, between 1995 and 2012. Shown is 
the complete set of maser transitions identified in the fully 
surveyed frequency interval between 18 and 26 GHz (with 
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1 Based on observations with the 100-m telescope of 
the MPIfR (Max-Planck-Institut fur Radioastronomie) at 
Effelsberg. 



Fig. 1. NH3 spectra showing a new maser component near 
a Local Standard of Rest (LSR) velocity of 45kms~ 1 . Note 
the slight velocity shifts in the (6,6) and (7,6) emission as 
a function of time (see Sect. 4.5). Less prominent quasi- 
thermal emission is centered near 60kms _1 . All spectra 
were taken at a B i 9 50 = 19 h 21 m 22 s .19, (5 B i950 = +14° 25' 
17'.% corresponding to a J20 oo = 19 h 23 m 39 S S3, S J20 qo = 
+14° 31' lO'.'l. For telescope beam sizes and channel spac- 
ings, see Table 1. 



the two orthogonal polarizations being averaged) and by 
observations of the (9,9), (10,10), and (12,12) transitions 
between 27.5 and 31.4 GHz. Lines, observing dates, channel 
spacings, and excitation above the ground state are sum- 
marized in Table 1. The characteristic rms noise level in a 
0.7-1.0kms _1 channel is of order 0.015 Jy. Note that the 
(3,1) line at 22.23450 GHz is blended by the much stronger 
61,6 — ► 52,3 line of H 2 0. 

How can we be sure that each of the spectral line- 
shapes displayed is affected by population inversion? A 
strong argument favoring maser emission is a comparison 
of all detected lines. While Figs. HHZl display data from a 
total of 19 NH3 inversion lines, a significantly larger num- 
ber of such transitions was actually detected in the rich 
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Table 1. Summary of NH3 maser observations 



Line 


V 


Epoch 


Channel 


-Eiow/k 


S 


Vlsr 


AV 1/2 










spacing 










(J,K) 


(GHz) 






(kms -1 ) 


(K) 


(Jy) 


kms 


- 1 


(6,2) 


18.884695 


2012, 


Apr.4-8 


0.97 


578 


0.22 


54.39±0.02 


1.26 


(5,3) 


21.285275 


2012, 


Apr.4-8 


0.86 


380 


0.30 


57.01±0.02 


1.37 


(6,3) 


19.757538 


2012, 


Apr.4-8 


0.93 


551 


0.35 


56.50±0.01 


1.44 


(5,4) 


22.653022 


2012, 


Apr.4-8 


0.81 


342 


0.52 


57.00±0.02 


1.54 


(7,4) 


19.218465 


2012, 


Apr.4-8 


0.95 


713 


0.09 


54.74±0.05 


1.31 


(7,5) 


20.804830 


2012, 


Apr.4-8 


0.88 


664 


3.10 


54.44±0.02 


0.95 


(8,5) 


18.808507 


2012, 


Apr.4-8 


0.97 


892 


0.14 


54.38±0.02 


1.44 


(6,6) 


25.056025 


2008, 


Oct. 25 


0.73 


406 


6.20 


45.72±0.01 


1.60 






2012, 


Apr.4-8 


0.73 




5.30 


47.22±0.01 


1.53 


(7,6) 


22.924940 


2008, 


Oct. 25 


0.80 


605 


1.95 


46.02±0.01 


1.77 






2012, 


Apr.4-8 


0.80 


605 


2.00 


47.55±0.02 


1.67 


(8,6) 


20.719221 


2012, 


Apr.4-8 


0.88 


833 


1.29 


45.92±0.06 


1.12 


(9,6) 


18.499390 


2012, 


Apr. 5 


0.99 


1089 


0.80 


52.25±0.45 


1.00 














6.80 


54.37±0.04 


1.19 














3.75 


56.80±0.04 


1.53 














1.15 


61.23±0.20 


3.70±0.50 






2012, 


Apr. 7 


0.99 




7.80 


54.28±0.03 


1.00 














4.90 


56.59±0.06 


1.02 














1.10 


60.97±0.24 


4.01±0.54 


(7,7) 


25.715182 


2008, 


Oct. 23 


0.71 


535 


1.70 


45.60±0.02 


1.46 








0.02 




2.48 


46.08±0.01 


1.27±0.03 






2008, 


Oct. 25 


0.71 




2.27 


45.66±0.05 


1.57 






2012, 


Apr.4-8 


0.71 




0.20 


47.29±0.11 


1.38 














0.18 


49.92±0.11 


1.23 


(9,7) 


20.735452 


2012, 


Apr.4-8 


0.88 


1020 


0.08 


54.72±0.08 


0.88 


(10,7) 


18.285434 


2012, 


Apr.4-8 


1.00 


1303 


0.03 


54.34±0.05 


1.00 


(9,8) 


23.657471 


2012, 


Apr.4-8 


0.77 


940 


0.22 


54.72±0.08 


0.93 


(9,9) 


27.477943 


1995, 


Nov. 28 


0.53 


848 


0.16 


43.76±0.07 


2.14±0.15 








0.13 




0.15 


43.79±0.11 


2.92±0.32 






1995, 


Dec. 19 


0.53 




0.12 


43.63±0.19 


2.70±0.41 






2011, 


Aug. 31 


0.94 




0.07 


46.81±0.04 


1.16 


(10,9) 


24.205287 


2012, 


Apr.4-8 


0.76 


1133 


0.34 


46.06±0.04 


1.20 


(11,9) 


21.070739 


2012, 


Apr.4-8 


0.87 


1446 


1.24 


47.12±0.01 


1.51 


(12,12) 


31.424943 


1995, 


Nov. 28 


0.46 


1452 


0.14 


43.84±0.24 


1.72±0.38 








0.12 




0.12 


44.30±0.30 


2.35±0.71 






1995, 


Dec. 19 


0.47 




0.10 


44.39±0.36 


1.77±0.79 






2011, 


Aug. 31 


0.06 




0.60 


46.87±0.02 


0.81±0.03 



The data are ordered with respect to K from 2 to 12 and, within a K ladder, with rising J (see Sect. 1). Column 5: Ei ovr /k: 
Energy above the ground state of the lower level of a given inversion doublet; k is the Boltzmann constant; (E up - i?i ow )/k ~ 
1.0-1.5 K. The beam size is 8b ~ 50" x (18GHz/^) with v being the line frequency. Column 6: Uncertainties in the flux densities 
are dominated by calibration errors, ± 15% (1<t; see Sect. 2) and are thus higher than the formal errors derived from Gaussian fits. 
Furthermore, if lines are spectrally unresolved, the given peak flux density is only a lower limit. Column 7: From Gaussian fits 
with formal standard deviations, which are likely underestimating the real errors in case of spectrally unresolved features. Column 
8: In many cases, the line width of the maser component is similar or smaller than the channel spacing displayed in column 4. In 
these cases, obtained line widths are upper limits, peak flux densities are lower limits, and no error is given. No maser emission 
was seen in the (1,1), (2,1), (4,1), (5,1), (6,1), (2,2), (3,2), (4,2), (5,2), (3,3), (4,3), (7,3), (4,4), (6,4), (5,5), (6,5), (8,7), (10,8), 
(12,9) and (10,10)) lines down to 3a limits of order 0.05 Jy in an individual 0.7-1.0 kms 1 channel. 



molecular spectrum toward W51-IRS2 (see the caption to 
Table 1). All these spectra, e.g. the metastable (J=K) 
(1,1), (2,2), ... (5,5), and (10,10) transitions exhibit the well 
known comparatively broad Vlsr ~ GOkms^ 1 feature (e.g., 
Mauersberger et al. 1987). Since its intensity is closely re- 
lated to excitation of the respective levels above the ground 
state (column 5 of Table 1 for the lines discussed in the fol- 
lowing) and since the feature is not time variable within 
the limits of accuracy over the last quarter of a century, 
the emission must be quasi-thermal. 

Each of the 19 NH3 transitions discussed in this paper 
shows in addition to this quasi-thermal component at least 



one additional spectral feature. These features tend to be 
narrower and almost all of those observed repeatedly show 
signs of variability (see Sects. 4.2 and 4.3). Furthermore, 
none of these additional components is seen in the inver- 
sion lines with lowest excitation (e.g., (1,1), (2,2), (3,3) or 
(5,5)), which were part of the frequency survey in 2012. All 
this strongly hints at maser emission, since such features, if 
quasi-thermal, would have to be seen in all inversion lines 
below a given excitation level. We thus adopt the maser 
classification for all 19 transitions displayed in Figs. [THTJ 
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Fig. 2. More NH3 spectra showing the new 45kms~ 1 com- 
ponent (see also the caption to Fig. [TJ . In the second low- 
est panel, the features near 30 and QOkms" 1 are quasi- 
thermally excited satellite lines of the (7,7) transition, 
caused by hyperfine splitting. 



3.2. The spectra 

Among the 19 NH3 maser lines identified by us in W51- 
IRS2, 13 have not been reported before from this source. 
The (6,3) and (9,6) lines (Madden et al. 1986), the (5,4), 
(7,5), and (9,8) lines (Mauersberger et al. 1987), and the 
(11,9) line (Wilson & Henkel 1988) were known before this 
study. However, the non-thermal features in the (6,2), (5,3), 
(7,4), (8,5), (6,6), (7,6), (8,6), (7,7), (9,7), (10,7), (9,9), 
(10,9), and (12,12) transitions are new discoveries in W51- 
IRS2, with the (12,12) transition being the most highly 
excited NH3 maser line so far detected. Among the newly 
identified maser lines in W51-IRS2 only the (6,6) and (8,6) 
transitions have an analog in NGC 6334 (Beuther et al. 
2007; Walsh et al. 2007) and possibly in IC 342 (Lebron 
et al. 2011). All other lines have, to our knowledge, not 
yet been reported to show maser emission in outer space. 



Fig. 3. Additional NH3 spectra showing the new 45kms 1 
component (see also Figs. [TJ and [5]). 



Figs. HHZl thus reveal eleven new molecular maser transi- 
tions. 

Figures [JJ [3] display spectra, which contain a previously 
absent (e.g., Mauersberger et al. 1987) NH 3 velocity com- 
ponent in W51-IRS2, near a Local Standard of Rest veloc- 
ity of Vlsr = 45kms~ 1 . In most cases the broad quasi- 
thermal feature near 60kms _1 is dwarfed, sometimes to 
such an extent, that it is only barely visible (e.g., Fig. [TJ). 
Figure 2] shows a five position map (spacing: 30") of the 
(J,K) — (7,7) maser. Averaging the offset positions, the 
amplitude of the resulting signal relative to that from the 
central position is 0.18±0.04. Using a Gaussian beam size 
of 0b = 35", we obtain for a point source, in good agree- 
ment, a corresponding ratio of 0.14. Observations of the 
nearby source W51e in the (7,7) line on October 23, 2008, 
showed no 45kms~ 1 maser component. The same holds 
for a (6,6) spectrum obtained on October 25, 2008, toward 
W51e. The 5a limits are 0.27 Jy and 0.88 Jy in 0.71 km s" 1 
and 0.73kms _1 channels, respectively. 
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Fig. 4. A map of NH3 (J,K) = (7,7) emission with a spac- 
ing of 30" toward W51-IRS2. The velocity scale is Local 
Standard of Rest (LSR), the ordinate gives flux density 
(Jy). North is up and east is left. For beam size and chan- 
nel spacing, see Table 1. 

Figure [5] displays spectra of the 45kms~ 1 component 
obtained with high spectral resolution (see also Sect. 4.5). 
While the (12,12) spectrum from November 1995 has been 
smoothed, the original data, albeit being noisy, indicate 
that the line is as broad as the (9,9) feature measured at 
the same epoch. Figs. [6] and [7] show spectra of ammonia 
lines with maser features being closer to the quasi-thermal 
"systemic" emission peaking near 60kms~ 1 . Interestingly, 
the (3,3) and (10,8) lines, previously reported to exhibit 
maser emission at ^54kms~ 1 (Mauersberger et al. 1987; 
Zhang & Ho 1995), now show exclusively the quasi-thermal 
60kms _1 component (Fig. [5]). While this requires a real 
change in the overall (10,8) line shape, the (3,3) maser was 
discovered through interferometric measurements and may 
be too weak and too close in frequency to the peak of the 
quasi-thermal emission to be detectable with single-dish 
telescopes. 

A comparison of the two linearly polarized components 
of the masers measured between 18 and 26 GHz shows no 
significant differences in line shape in almost all cases. The 
exception is the (9,6) line, where the 57kms~ 1 component 
appears to be highly polarized (see Fig. [9]). An alterna- 
tive explanation, that there exist two maser sources and 
that there is a beam squint between the two orthogonal 
polarisations, causing the different line shapes, can be ex- 
cluded. Pointing and calibration measurements (Sect. 2) do 
not show any such position offset. 
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Fig. 5. High velocity resolution spectra of NH3 maser tran- 
sitions showing the new ^45kms _1 component. The upper 
and central profiles are extracted from Figs. [5] and 02 re- 
spectively. The second panel from the top shows a profile 
with eight contiguous channels smoothed (channel spac- 
ing: 0.96kms _1 , otherwise, see Table 1). The original spec- 
trum, albeit noisy, indicates a line width consistent with 
that shown in the spectrum of the (9,9) line measured on 
the same day. The Gaussian fit parameters of the (12,12) 
spectrum from November 28, 1995, given in Table 1, were 
deduced from the original unsmoothed spectrum. 



4. Discussion 

4.1. General considerations 

We have covered the entire 18-26 GHz frequency band, 
which includes all ammonia inversion doublets with low ex- 
citation. Nevertheless, when employing a 9b ~ 40" beam, 
Table 1 unambiguously tells us that only highly excited am- 
monia inversion doublets exhibit velocity components with 
inverted populations. Among the 19 detected maser lines, 
the (J,K) = (5,4) line at E ~ 340 K is closest to the ground 
state. 



Another remarkable feature is the number of ortho- ver- 
sus para-NH3 maser lines. Even though para-NH3 inversion 
doublets (K = 1, 2, 4, 5, 7, 8...) are twice as numerous as 
ortho-NH 3 doublets (K = 3, 6, 9...; the K = ladder is 
characterized by singlets), we get ten ortho- and only 9 
para-NH3 maser transitions. This may be understood in 
view of two important differences between the ortho- and 
para-species. Firstly, the ortho levels have twice the statisti- 
cal weight, thus providing larger column densities and pos- 
sibly higher maser amplification factors for ortho-to-para 
abundance ratios of order unity (e.g., Umemoto et al. 1999; 
Goddi et al. 2011). Secondly, as already mentioned, the 
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Fig. 6. NH3 maser spectra taken from W51-IRS2 not show- 
ing the 45kms _1 maser component. For beam sizes and 
channel spacings, see Table 1. 



K=0 ladder contains single rotational states. In the gas 
phase, ortho- and para-Nti3 can exchange only very slowly 
(e.g., Cheung et al. 1969). The ortho-NH 3 K = single 
states therefore connect to one of the two levels in a K = 3 
inversion doublet, thus providing potentially significant de- 
viations from quasi-thermally balanced populations at K — 
3 (e.g., Walmsley & Ungerechts 1983). Within this context 
it is worth mentioning that the two strongest maser lines, 
those from the (6,6) and (9,6) inversion doublets, belong to 
the ortho-species. 

Another noteworthy fact is that among the 19 maser 
lines detected by us, only four are metastable (J—K): three 
ortho-NH 3 [(6,6), (9,9), and (12,12)] transitions and one 
para-NH 3 [(7,7)] line (see Table 1). The vast majority of 
the maser lines detected in this study originates from non- 
metastable inversion doublets (J>K). Their populations 
decay rapidly to the metastable levels, so their excitation 
is complex and either requires extremely high densities 
(<;10 6 cm~ 3 ) and/or an intense infrared radiation field. 
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Fig. 7. NH3 maser spectra taken from W51-IRS2 not show- 
ing the 45kms~ 1 maser component (see also Fig. 



4.2. Variability on short time scales 

We seached for changes in the line profiles of our spectral 
broadband survey during April 4-8, 2012. The (5,3), (5,4), 
(7,5), (7,6), (8,6), (9,7), and (11,9) transitions remained 
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stable, within the calibration uncertainties (Sect. 2), over a 
time interval of four days. For the (6,3), (7,4), (6,6), and 
(7,7) lines the same holds for an interval of three days, while 
we see no significant change in the (6,2), (8,5), (9,8), and 
(10,9) transitions (and in the (7,7) transition in October 
2008, see Figs. [5] and H|) after two days. However, the (9,6) 
line shows a significant variation in line shape within a time 
interval of only two days (Fig. [7]) . To our knowledge this is 
the fastest variation so far seen in an ammonia maser and 
might indicate that the maser is due to exponential am- 
plification of radiation, that is unsaturated. Since it is one 
of the two strongest ammonia masers toward W51-IRS2, 
we could then speculate that all the 19 masers we have 
found are unsaturated. However, data with higher angular 
resolution are definitely needed to further clarify this point 
(see, e.g., the arguments of Wilson et al. 1991 favoring a 
saturated (9,8) maser line). 

We can assume that the time scale for intensity varia- 
tions is equal to the light travel time across the source. The 
main variation in the (9,6) line is seen near 55kms~ 1 . If 
we take a flux density of S ~ 2 Jy (Fig. O, a beam size of 
48 '.'6 (see caption to Table 1), a flux to main beam bright- 
ness temperature conversion factor of 1.5, and a distance 
of 5kpc (see Sect. 1), we then obtain an angular scale of 
order O'/l and a corresponding brightness temperature of 
order 10 6 K, similar to the limit obtained for the (9,8) line 
by Gaume et al. (1993). 

4.3. Variability on long time scales 

Beside the Vlsr ~ 62kms~ 1 component of the (J,K) = 

(9.6) transition (Fig. [9]), which is too intense and time 
variable to be of quasi-thermal origin (cf. Madden et al. 
1985; Wilson et al. 1988), there are three additional fam- 
ilies of maser lines. These show emission near 57kms _1 , 
near 54.5km.s~ 1 , and near 45kms _1 (see Table 1). In the 
following we discuss all components in order of decreasing 
velocity. 

High velocity features (<;60kms ) were observed by 
Madden et al. (1986) and Wilson & Henkel (1988) in W51- 
IRS2. The Madden et al. (9,6) spectrum from late 1984 or 
early 1985 shows about half a dozen of narrow spikes be- 
tween 50 and 62kms _1 . Flux densities are a few Jy. In this 
case, however, the angular resolution is 96", so that the 
neighboring W51e region is also inside the beam (see, e.g., 
Fig. 1 of Gaume et al. 1993). In September 1987, Wilson & 
Henkel (1988), using the Effelsberg telescope and therefore 
exclusively measuring the W51-IRS2 region, find a partic- 
ularly strong component (^10 Jy) near Vlsr ~ 65kms , 
which is, however, absent in July 1988. 

Vlsr ~ 57kms _1 masers are seen in the (5,3), (6,3), 
(5,4), and (9,6) transitions, i.e. mostly in ortho-NH3. There 
is rapid variability in the (9,6) line (see previous subsec- 
tion). In the eighties (Madden et al. 1986; Wilson & Henkel 
1988), (6,3) emission was seen at 54kms~ 1 . The line has 
switched from one group of maser components to another. 
The (5,4) feature has varied little during the past quarter 
of a century, while maser emission from the (11,9) line, de- 
tected between August 1986 and October 1986 (Wilson & 
Henkel 1988), is not present in more recent spectra. 

The 54.5kms _1 velocity component is found in eight 
maser lines, namely in the (6,2), (7,4), (7,5), (8,5), (9,6), 

(9.7) , (10,7), and (9,8) transitions. In earlier studies, also 
(3,3), (10,8), and (11,9) maser emission has been reported 
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Fig. 8. The two previously reported maser lines from W51- 
IRS2 (Mauersberger et al. 1987; Zhang & Ho 1995), which 
were found to only exhibit quasi-thermal emission. The 
triple component line shape of the NH3 (3,3) transition 
is caused by hyperfine splitting. Channel spacings are 
0.88 kms" 1 for the (10,8) and 0.77kms- 1 for the (3,3) pro- 
files, respectively. 



(Mauersberger et al. 1987; Wilson & Henkel 1988; Zhang 
& Ho 1995). The (6,3) line, which 25 years ago also con- 
tributed to the 54.5 km s -1 group of masers (Madden et al. 
1986), is now seen at ~57kms _1 (see above). Except the 

(9.6) transition, all lines belong to the para-NH3 species. 
Mauersberger et al. (1987) searched for (6,2) emission but 
the line remained undetected. Their upper limit, 0.1 Jy, is 
well below the 0.3 Jy signal we have found. The (7,5) and 

(9.7) transitions did also not show a maser component (see 
Fig. 1 of Mauersberger et al. 1987). However, in July 1988, 
Wilson & Henkel (1988) detected a 0.4 Jy maser feature 
in the (7,5) line. Since then, the maser peak flux density 
has increased sixfold. The opposite trend, decreasing in- 
tensities with time, is seen in the (11,9), (9,8), and (10,8) 
lines. (11,9) maser emission was observed in 1986 (Wilson 
& Henkel 1988) but has faded. The (9,8) and (10,8) lines 
were much stronger ^25 years ago. For the (9,8) line, al- 
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Velocity (km/s) 

Fig. 9. Linear polarisation in the 57km s component of 
the (9,6) maser line on April 5 and 7, 2012. Shown are 
the two orthogonal components. The left panels display the 
horizontal (approximately east-west) and the right panels 
the vertical (approximately north-south) component, be- 
cause the line was observed close to transit. Since the line 
is strong, the difference is highly significant and is visible 
in each of the 20 (April 5) and 60 (April 7) individual sub- 
scans. For beam size and channel spacing, see Table 1. 



ready Wilson et al. (1991) reported a moderate decrease by 
^30% relative to measurements during the previous years. 
Now the decline has reached 90% and the previously ob- 
served (10,8) maser is now below our detection threshold. 

Another eight transitions, the (6,6), (7,6), (8,6), (7,7), 
(9,9), (10,9), (11,9), and (12,12) lines, belong to the V LS r ~ 
45kms~ 1 group. Almost all of these are ortho-NH3 transi- 
tions, again with one notable exception, this time the (7,7) 
line. The 45kms _1 component is the only one not seen 
in the (9,6) maser transition. For the broad quasi-thermal 
(8,6) feature near 60kms _1 Mauersberger et al. (1987) pro- 
posed anti-inversion, implying a low but still positive ex- 
citation temperature and a high optical depth. It is thus 
remarkable that we find now, near 45kms~ 1 , strong (8,6) 
maser emission. In the (9,9) and (12,12) lines, which have 




(J.K) 



Fig. 10. Selection rules for purely rotational transitions be- 
tween two NH3 inversion doublets belonging to the same 
K-ladder. Parities are given on the right hand side. 



been repeatedly measured, there is a tendency of rising peak 
flux densities in more recent observations. On the contrary, 
in spring 2012, the (7,7) line has dramatically weakened, re- 
vealing a previously unseen second maser feature (Fig. [5]). 

4.4. Constraints for pumping scenarios 

We noted in Sect. 4.2 that all observed masers in W51-IRS2 
might be unsaturated. Wilson et al. (1991) proposed in- 
stead saturated emission in the case of the (9,8) line. The 
two arguments were a large line width (~0.4 km s" 1 ; Wilson 
& Henkel 1988) and constant peak flux density. In the 
meantime, however, the line has varied substantially and 
the line width is rather narrow. Measuring the line widths 
of all detected masers and comparing these, if possible, with 
previously taken data with high spectral resolution will re- 
veal the correlation between line width and peak flux den- 
sity and will thus provide a definite answer to the question 
of maser saturation. This is part of a follow-up study. 

Figure [TU] illustrates an important rule for radiative 
dipole transitions within a K-ladder containing NH3 inver- 
sion doublets. The upper level of the upper doublet (the 
(J+l,K) inversion doublet) is directly connected to the 
lower level of the lower doublet {J,K) and vice versa. Thus, 
if the upper doublet is inverted, the lower doublet should 
be anti-inverted. 

This is not confirmed by our measurements. The maser 
emission at 57kms _1 is observed in the (5,3) and (6,3) 
transitions. The 54.5 km s -1 component is seen in the (7,5) 
and (8,5) and in the (9,7) and (10,7) lines. In case of the 
45kms~ 1 component, we even find the (6,6), (7,6), (8,6) 
and the (9,9), (10,9), (11,9) transitions. It appears as if, 
contradicting the above mentioned rule, the combination of 
(J+l,K) and (J,K) masers is not forbidden but actually 
preferred. The former presence of the (9,8) and (10,8) maser 
lines and their simultaneous weakening also points into this 
direction, even though truly systematic variability studies 
of such maser pairs are still not available. 

As already noted by Madden et al. (1986), there are 
three main pumping schemes which could, in principle, ex- 
plain the observed NH3 maser emission. These are (1) line 
overlap, (2) pumping by infrared radiation from the dust, 
and (3) pumping by collisions. Simulations of NH3 maser 
emission encounter several problems. One is that collision 



8 



Henkel, C, Wilson, T.L., Asiri, H., Mauersberger, R.: New ammonia masers 



3 ^(J+ 1.KJ 
+ 

(J,K) 

Fig. 11. Pseudo-selection rules for rotational transitions 
between two NH3 inversion doublets belonging to the same 
K-ladder, when vibrational excitation is involved. Parities 
are given on the right hand side. Excitation into a vibra- 
tionally excited level and subsequent decay leads to A J = 
0, ±1, and ±2 pseudo-transitions without parity change in 
the ground state. Shown are two R-branch transitions with 
AJ = -1. 



rates are only known for inversion levels J < 6 (e.g., Danby 
et al. 1988). Another is the spectrum of the infrared radia- 
tion field. For vibrational excitation of NH3 by IR photons 
near A = 10^m (see Bally et al. 1987), the silicate dust fea- 
ture may significantly perturb the radiation field, making 
any prediction highly speculative (Brown & Cragg 1991). 
Madden et al. (1986) suggested for the NH 3 (10,6) + -(9,6)~ 
and (10,6)~-(9,6) + transitions, the upper indices denoting 
parity, line overlaps. Such blends, if they are not common 
and only affect the (9,6) line, may explain the very par- 
ticular behavior of this inversion line. Each maser family 
(Sect. 4.3) is mostly composed of either ortho- or para-NH3 
transitions. This indicates that the two ammonia species 
are, as expected, reasonably well separated. 

Schilkc (1989) extrapolated NH3 collision rates to high- 
J levels and even included vibrational excitation (the 1*2 = 
1 level) , finding that vibrationally excited levels strongly af- 
fect the excitation of ground state inversion lines. For T^ust 
~ Tkin, excitation temperatures can raise well above T^ust 
and Tkin, in particular when temperatures above 160-180 K 
are being reached. The latter does not hold for the bulk of 
the NH3 bearing gas in Orion-KL (Hermsen et al. 1988) 
but is realized in W51-IRS2 (Mauersberger et al. 1987). 
Interestingly, Schilke (1989) mentions the (9,8) line, which 
is easily inverted under such conditions, while the (10,8) 
line should not become a maser. This is consistent with our 
data from April 2012 (Figs. [7] and H]), but inconsistent with 
previous data also showing the (10,8) line to be inverted 
(Mauersberger et al. 1987). 

Brown & Cragg (1991) also studied the influence of vi- 
brational excitation on the ground state transitions. They 
could reproduce (6,3) maser emission and predicted quasi- 
thermal emission for the (4,3) and (5,3) lines. While the 
former is consistent with our data, the latter is not (see 
Fig. ©, so that our data do not fully confirm their model 
in the case of W51-IRS2. 

Following Schilke (1989), vibrational excitation leads to 
pseudo-selection rules in the ground state, where then A J 
= 0, ±1, and ±2 "transitions" are allowed without parity 



change. This directly contradicts the selection rules illus- 
trated in Fig. [TOl and leads to the situation in Fig.[TTJ which 
predicts maser emission in several adjacent J levels within 
a K-ladder, as observed. We conclude that the measured se- 
ries of masers in (J,K) and ( J+1,K) lines, sometimes also 
including ( J+2,K), are a direct consequence of vibrational 
excitation. Such a scenario was already proposed for the 
(3,3) and (4,3) lines of 15 NH 3 in NGC 7538-IRS1 (Schilke 
et al. 1991). However, our case with several pairs and even 
triples of maser lines within a given K-ladder is a much 
more convincing one. The apparent contradiction to Schilke 
with respect to the (10,8) line may be a result of insuffi- 
cient knowledge of collisional parameters (available only up 
to the (6,6) line; Danby et al. 1988), spatial fine structure, 
and the radiation field. 



4.5. A//-/3 masers, the ambient molecular cloud, and a feature 
with velocity drift 

Interferometric measurements to determine accurate posi- 
tions have been carried out for the (3,3) and (9,8) maser 
lines (Gaume et al. 1993; Zhang & Ho 1995). The (3,3) and 
(9,8) positions (a J200 o ~ 19 h 23 m 39 '.'8, <5 J20 oo ~ +14° 31' 
05 '.'0) agree within the uncertainties of up to one arcsec- 
ond and may represent the position for the entire family of 
54.5 km s -1 maser lines (see Sect. 4.3). These are associated 
with an ultracompact HII region (Wilson et al. 1991) that 
is also known to show 6.7 GHz methanol emission (Phillips 
& van Langevelde 2005). For the family of 57kms _1 masers 
an accurate position is not yet known. 

While frequently seen in late-type stars, W51-IRS2 is 
one of only three star forming regions known to show 
SiO maser emission (Snyder & Buhl 1974; Hasegawa et al. 
1986). Such emission originates from vibrationally excited 
levels more than 1000 K above the ground state and thus 
requires highly excited gas (e.g., Goddi et al. 2009). There 
are no high NH3 resolution data to prove a spatial correla- 
tion between ~45kms~ 1 NH3 and SiO emission. However, 
there is a 45kms _1 component seen in the SiO v—2 J = 
1— >0 maser line, possibly tracing an accelerating bipolar 
outflow close to the obscured protostar W51-North (Eisner 
et al. 2002; Zapata et al. 2009). The possible correlation 
between the highly excited 45 km s _1 v — 2 SiO maser com- 
ponent (a J20 oo = 19 h 23 m 40 '.'06, 6 J2 aoo = +14° 31' 05 "6) 
and NH3 is supported by the fact that the ^45kms _1 
NH3 masers are those with highest excitation, including 
the (12,12) transition at E/k ~ 1450K (Table 1). 

Also related to this component, and perhaps the most 
striking aspect of our observations, is the velocity drift of 
the measured ^45kms _1 NH3 features with time. This is 
seen in Figs. HH3] and even more so in Fig.[S] The drift is of 
order +0.2kms~ 1 yr^ 1 (Fig. [12]). For July 1994, extrapo- 
lating the dashed line connecting data from 1995 and 2012, 
we find a velocity of Vlsr ~ 43.5 km s -1 . For July 1998, 
Vlsk ~ 44.5 km s . Both values correspond to distinct SiO 
maser velocity components (Eisner et al. 2002), detected at 
these two epochs and might be related to the lower velocity 
SiO feature monitored by Fuente et al. (1989; their Fig. 2), 
which shows a redward drift of order 0.4-0.5 km s -1 yr -1 
between 1985 and 1989. In 1994, the NH 3 maser velocity 
is consistent with the strongest SiO maser feature observed 
by Eisner et al. (2002), if we neglect the potential influence 
of missing flux in these high resolution data. All this sup- 
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Fig. 12. Radial velocity versus epoch for the Vlsr ~ 
45 kms -1 component with the data from Table 1. The 
dashed line provides an approximate drift, based on the 
average velocity obtained in November 1995 and in April 
2012. For the latter epoch, the high velocity feature (Vlsr 
= 49.9kms _1 ) in the (7,7) line (Fig. [5]) was considered as 
an outlier and was not part of the average. 



The same holds for the circumstellar envelopes of late- type 
stars, which do not show NH3 maser emission. This implies 
that we have detected a new phenomenon which might be 
rare. However, the number of known star forming regions 
exhibiting SiO maser emission is still too small to provide 
a sound statistical basis. 



4.6. Why does W51-IRS2 host so many ammonia maser 
sources? 

A question still to be addressed is why does W51-IRS2 
host so many NH3 masers? One potentially useful property 
may be a fortuitous spatial configuration, with the dens- 
est part of the molecular hot core being located in front of 
an H11 region (this may, however, not be the case for the 
45kms _1 component if it is associated with the SiO source 
(Eisner et al. 2002)). Another ingredient is likely the ex- 
tremely high column density of NH3 (Sect. 1), which cer- 
tainly facilitates the build up of sufficient maser intensity. 
This holds in particular in the case of unsaturated maser 
emission (see Sects. 4.2 and 4.4), where larger column den- 
sities lead to exponential effects. Possibly most important, 
however, is the kinetic temperature of the bulk of the am- 
monia emitting gas. With T^ in ~ 300 K (Mauersberger et 
al. 1987), it is much higher than, for example, in Orion-KL 
(Hcrmsen et al. 1988). This high temperature at densities, 
where dust and gas temperatures are likely coupled, facili- 
tates vibrational excitation by infrared photons near 10 fim 
(see Sect. 4.4). This may cause significant deviations from 
quasi-thermal conditions in the ground state (Schilke 1989; 
see also Mauersberger et al. 1988a and Schilke et al. 1990, 
1992 for vibrationally excited NH3 in the cooler Orion-KL 
hot core). Line overlaps may also play a role. Note that the 
excitation above the ground state of the NH3 (12,12) maser 
line (Figs. [2] and [5]) is of the same order of magnitude as 
that of the first vibrationally excited state. 



ports a spatial correlation between SiO and the 45 km s 
NH3 maser component. 

From SiO, Eisner et al. (2002) have identified a linear 
structure with a northern (more blue-shifted) and southern 
(more red-shifted) arm along a position angle of ^105° with 
a velocity gradient of about 0.1 kms" 1 AU _1 along each 
arm. Assuming that the 45 kms -1 NH3 feature is an out- 
flow region or a bullet belonging to this structure, the NH 3 
velocity drift would correspond to a motion of 2 AUyr -1 , 
roughly eastwards along P. A. ~ 105°. 

Eisner et al. (2002) interpreted the SiO distribution in 
terms of an accelerating outflow from a deeply obscured 
massive stellar source with a 3-dimensional acceleration 
of 0.5 kms -1 yr _1 and 0.2kms _1 yr _1 projected along the 
line-of-sight. The latter agrees remarkably well with our 
ammonia data (Fig. IT21 . The agreement is so good that 
we refrain from suggesting any other model as long as fu- 
ture interferometric NH3 measurements are not indicating 
a clearly different scenario. 

Considering the proposed association of SiO and NH3 
masers, we may ask whether such a phenomenon may be 
common. Toward W51-IRS2, one out of several NH3 maser 
velocity components appears to be related. Toward the 
other two star forming SiO maser sources, Orion-KL and 
Sgr B2 (Hermsen et al. 1988; Huttemeister et al. 1993; 
Goddi et al. 2011), such NH3 masers are not observed. 



5. Conclusions 

Using the 100-m telescope at Effelsberg, we have searched 
for NH3 maser emission in the prominent massive star form- 
ing region W51-IRS2. The main results from this study, 
covering 16 years of monitoring, are briefly mentioned: 

— A total of 19 NH3 maser lines was detected. Eleven 
of these inversion transitions have been identified as 
masers for the first time in the interstellar medium. 

— The (9,6) line has been found to vary significantly 
during a time interval of only two days. This is the 
fastest change so far detected in an ammonia maser. 
Furthermore, its 57kms~ 1 velocity component shows 
strong linear polarisation, which is not found in any 
of the other lines. While variations in other transitions 
are slower, not only ortho- but also para-NH3 transi- 
tions show drastic changes during the last 16 years. The 
(12,12) line, at ^1450 K above the ground state, is the 
most highly excited ammonia maser line so far detected. 

— Apart from the (9,6) line, which is unique among the 
19 transitions reported here, there are three families of 
maser transitions, those exhibiting emission at 57, 54.5 
and 45 kms" 1 . The latter family was not known previ- 
ously. 

— There are many maser pairs in the same K ladder, such 
as (J,K) and (J+1,K). In some cases, this extends to a 
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triple, including (J+2,K). This observation gives strong 
support for vibrational excitation, induced by A ~ lO^tm 
photons, as the main pumping mechanism. 

— The 45kms~ 1 component masers show a secular drift 
of about 0.2 km s" 1 yr -1 and may arise close to the SiO 
maser source. Possible scenarios explaining the drift in 
terms of outflowing gas or orbital motion are outlined. 

— The most likely reason that W51-IRS2 hosts many more 
NH3 masers than any other known source is the com- 
bination of a very high temperature (~300K, ideal for 
10/im photons needed to vibrationally excite NH3), an 
exceptionally high NH 3 column density, and possibly 
also a fortuitous location with respect to continuum 
sources such as hyper-compact Hn regions. 

Besides measuring accurate line widths of many of the 
newly discovered maser lines to discriminate between satu- 
rated and unsaturated maser emission, the most important 
follow-up will be high angular resolution measurements. 
Only these will allow us to determine precise positions, 
which are essential to interprete the data in a broader con- 
text, also including other astrophysical tracers, and to ex- 
plain the origin of the enigmatic velocity drift. 
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